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Abstract — In this paper, we propose a new path planning method by using an electromagnetism-like mechanism algorithm. 
We use the different encoding method to solve trade-off problem that is encountered in the traditional path planning method. 
By combining different encoding method with electromagnetism-like mechanism algorithm, path point can be generated 
without the trade-off problem. In order to connect these points into smooth cur\>e, we compare two path smoothing method, 
Bezier curve and cubic splines interpolation. Finally, cubic splines interpolation is used because it can generate smoother 
path. 
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I. Introduction 

For all of the mainstream robots which have the movement capability, path planning will be a necessary and worthwhile 
problem to solve. The traditional path planning method is proposed by computer scientist Dijkstra [1], then A* algorithm [2] 
solves the problem of direction. The A* algorithm makes it work more efficiently, but these algorithms will encounter the 
map trade-off problem in the practical application. In a known environment, the map needs to be cut into many grids, and 
these grids are used to plan the path. 

The process of grids is a trade-off problem. If the map is cut into large grids, the relative generated path will be very rough. 
However, if the map is cut into small girds, it will cost a lot of computing costs. In real case, the obstacles on the map is 
irregular shape, such that path planning result is not the shortest path. As a result, many algorithms were proposed to 
compute shorter or more effective path for the trade-off problem. 

To solve the problem mentioned above, we need to change the processing method for path planning by using a different map 
encoding method [3] and connect with the heuristic algorithm to solve the path planning problem. The electromagnetism-like 
(EM) algorithm [4] is innovative in the heuristic algorithm which is based on the characteristics of electromagnetic attraction 
and repulsion. 

In this paper, instead of applying EM algorithm for optimization, we combine EM algorithm with different encoding method 
for path planning. It can avoid traditional trade-off problem by spreading points and connecting these points. The algorithm 
computes the path with path smoothing method. In this paper, cubic spline interpolation is chosen because it can smooth the 
path without increasing the algorithm’s parameters. To demonstrate the proposed method, we simulate the possible map 
configuration, and use map simulation to verify the feasibility of the algorithm. 

II. Electromagnetism-Like Mechanism Algorithm 

Electromagnetism-like Mechanism algorithm [4] is a heuristic algorithm proposed in 2003, and the convergence is presented 
in [5], EM algorithm applies the concept of Coulomb's law in the electromagnetic, the better electronic attract bad electrons, 
and bad electronics will repel other good electronics. The attraction and repulsion forces between electrons accumulate into 
total force in order to make the path result to better direction. Because EM algorithm has good computing ability, there are 
many kinds of applications, such as optimization question for communication system [6] [7], optimization the neural 
networks [8], optimization the learning rule for fuzzy neural networks [9], path-tracking problem [10], optimized tool path 
planning [11], parameter optimization [12], the routing problem [13], optimization the knapsack problem [14] and other 
applications [15] [16]. 

Besides considering the better solution, EM algorithm also consider the worst solution as a contrarian indicator, these 
multiple reference can let the solution converge faster. On the other hand, the EM algorithm has fewer parameters; it is more 
intuitive and convenient in parameter adjustment. And it can add more electronics in the map, such that user can adjust 
parameter based on road condition. It is feasible to choose the EM algorithm for path planning. 
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The EM algorithm includes four parts [17]: initialize the algorithm (Initiation), search the neighborhood to avoid falling into 
the optimal solution (Local), calculate the total force (CalcF) for each electronics, and move the electronics along the force 
direction (Move). 

General framework is shown below. 

1: Initiation() 

2: ite<—T 

3: while ite <MAXITE do 
4: Local (LITE,a) 

5: F <—CalcF() 

6: Move(F) 

7: ite <—ite + 1 
8 : end while 


In this paper we use the modified EM algorithm which adds the rearrange function in original EM algorithm. Rearrange 
function can adjust the position of electrons without altering the path, and it makes the EM algorithm can generate global 
optimal solution quickly. The flow chart of the EM in this paper is shown in Fig. 1. 



^_ 

Fig. 1. Flow Chart in This Paper 


Stepl: Initiation 

Initiation spreads point randomly in the map. Each one of the point is regard as one electrical point, and all of the electronics 
are connected from start point to the end point, as shown in Fig. 2(a). Then the point is spread, and each electronic evaluates 
the other electronics’ magnetic which has same dimensions, as shown in Fig. 2(b). 
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Fig. 2. (a) Link the Electrical Point. 
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Step2: Local 

Local makes a provisional electronics for each dimension. The electronic is a replica of the original one, then move 
randomly. If Local finds the better result, it will replace the original one, and finally find out whether there is a better global 
optimal solution. 

Step3: CalcF 


It is necessary to calculate charge amount for each electrons <2, before calculating the total force. The attraction and 
repulsion forces are determined by charge amount. The charge amount of the i th electron is shown in (1). 




(i) 


where m means the number of particles, e means the current electronic, and e besl means global optimal solution. The 
global optimal solution has a greater attraction. Then, calculate the total force F' for each electron, as shown in (2). 

QQ j 


F' = I. 
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Although equation (2) is more accord with physical meaning, it costs more time in algorithm. By [18], we use equation (3) to 
accelerate the implementation rate of algorithms. 

\(e i -e i )Q‘Q J if /(r'X./V) 

\fe‘-e ] )Q‘Q J if f(e J )>f(e') 

Step4: Move 


*"=Z 


(31 


After calculating the total force F ‘, electrons are designed to move a random distance based on the direction of the total 
force, and the movement method is shown on equation (4). 


4 =4 + YT—v i ( RNG '> 

\\r 


(4) 


where y means a uniform distribution between 0 and 1. It also can use other types of distributions in calculating the step 
length. RNG means a random number for step length. It can determine a nonzero probability to move to unvisited search 
space. 

Step5: Rearrange 


If the electronics are overly concentrated, the EM algorithm will occur sharp point phenomenon when it applied to path 
planning, as shown in Fig.3. Therefore, we need the Rearrange function to redistribute the electrons. As a result, the global 
optimal solution is generated more quickly. 



Fig. 3. Sharp Point Phenomenon 













International Journal of Engineering Research & Science (IJOER) 


ISSN, [2395-6992] 


[Vol-3, Issue-10, October- 2017] 


The general framework of the Rearrange function is described as follows. 

1: y <- RNG(QJ) 

2: if / < T * del then 

3: j<- 1 

4: for i = 1 to h do 

5: if length(i) = ——— 

n 

6: e[ = position(i) 

7: if i = n + 1 

8 : break 

9: end if 

10: j<~j +1 

11: end if 
12: end for 


y means a uniform distribution between 0 and 1. T means the threshold which makes the path rearrange, det is a parameter 
decrease when iteration number increases. We use det to control the probability which let algorithm execute the Rearrange 
function. L is the length of the path, n is the number of sections. 

III. Path Smoothing 

The result of the path planning often has a non-smooth turning characteristic. It is caused by the rough segmentation of the 
map. This makes the path planning result possess sharp angle turning in the actual application. Therefore, the vehicle must 
need a lot of deceleration, or even stop the vehicle to make a turn. 

After the path smoothing process, it makes the smooth path planning result. Therefore, the vehicle can make a turn without a 
lot of deceleration in actual application. 

However, when we use high-order polynomial to approximate the path planning result, Runge's phenomenon will occurs. 
Runge's phenomenon occurs when using polynomial interpolation with polynomials of high degree over a set of 
equispaced interpolation points. It is a problem of oscillation at the edges of an interval. To avoid Runge’s phenomenon, 
multi-sections are used rather than high-order interpolation method for path smoothing. This paper will introduce two 
methods of smoothing path in follow. 

3.1 Bezier Curve [19] 

Given point D 0 and D l , the Bezier Curve L 0 is a straight line between those two points. The curve is given by equation (5): 

L 0 =(l -t)D 0 +tD l (5) 

If Bezier Curve starts from D 0 to D , and has a relay point D ] , the second order Bezier Curve C 2 (t) is shown as equation (6): 

C 2 (r) = (1 - 1) 2 D 0 + 2t(l - t)D x +t 2 D 2 (6) 

The general form of a n degree Bezier curve defined by the control points /> (where i = 0, • • •, n ) is 


c( o=ZB,,„m 


where B Jt) is called Bernstein polynomials that are defined in (8). 






(7) 


( 8 ) 


the cubic Bezier curve C 3 (f) is given in (9). 

C 3 (t) = (1 - tfD 0 + 3rd - r) 2 D, + 3r (1 - t)D 2 + t 3 D 3 


(9) 
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3.2 Cubic Splines Interpolation [20] 

Spline interpolation uses segmented polynomial to interpolation, and it can achieve the smaller interpolation error by using 
low-order polynomial, and avoid Runge's phenomenon. 

First, define k r = C (t [ ) for i = 1,2,— 1, C (f ( .) = w j , C,(f i+1 ) = w i+1 and k 0 =0,k n = 0 . And it Lagrange form for C (t i ) is 


C'fy) k <+' (x ~ f ^ | 


Then the C, (x) and C t (x) are given as 


2«,2«„,-',) 

where A and II are shown in equation (13) and (14). 
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Therefore, C, and C, are generated as following equation (15) and (16). 
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Also, C (r) requires the following condition. 

c:_ 1 (f,) = c:(r,) i = 1,—1 

Then we can get equation (17). 
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By using =t i+1 —t i ,and d j = —^- L , equation (17) has 


C(t) = -—sk^, -—sk. +d. 
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By C,_!(f,) = C,(f,), then equation (20) is given. 
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Bezier curve is widely used in computer graphics and other path smoothing method. It is a simple and practical way. 
However, distinguishing between control point and path point in path planning algorithm is more complicated. 

In the EM algorithms, cubic splines interpolation is used in doing path smoothing algorithm so that the number of parameters 
does not increase. Regarding each point as one path point, the algorithm enhances the relevance and thus enhances the 
computing performance. 


IV. Simulation Result 

The simulation parameter is shown on Table. 1, and the simulation result is shown in Fig. 4. As shown in Fig. 4, the square is 
the start point, the star is the end point, and the circle is the obstacle. After the path planning and smoothing, the curve is the 
algorithm result. We calculate the cost function by using the Pythagorean Theorem. When the path contacts the obstacle, its 
cost function will be multiplied 100 as a penalty. Thus we can always find the best cost without touching the obstacle. 


Table 1 

Simulation Parameters 


Parameters 

Max 

Iteration (MAXITE) 

Local 

Iteration(/-/77i) 

Number of 
population 

Number 

of 

Control 

Point(m) 

Delta! det ) 

Rearrange 

threshold(T) 

Penalty 

Values 

50 

100 

10 

5 

0.4 

0.1 

100 


8 
7 
6 
S 
4 
3 
2 
1 
0 
-1 
-2 

-4 -2 0 2 4 6 8 

Fig.4. Simulational Result 
V. Conclusions 

This paper has presented a path planning algorithm to solve the trade-off problem that is encountered in the traditional path 
planning method. Comparing two path smoothing method, Bezier curve and cubic splines interpolation, cubic splines 
interpolation is selected to smooth the path without increasing the number of parameters. By using the path smoothing 
method, the vehicle could go through the angle smoothly like reality. 
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